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ABSTRACT: In this paper, we introduce an approach to increase the density of field-effect transistors 
in the framework of a wideband current amplifier with DC-offset cancellation. In the framework an 
approach we consider in the manufacturing of inverter in heterostructure with specific configuration. 
Several required areas in the heterostructure should be doped such as diffusion or ion implantation. 
After that the dopant and radiation defects should by annealed in framework optimized scheme. We 
also consider an approach to decrease the values of mismatch-induced stress in the considered 
heterostructure. We also introduce an analytical approach to analyze the mass and heat transport of 
heterostructures during the manufacturing of integrated circuits with account mismatch-induced 


stress. 
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INTRODUCTION 


In the present days several actual problems of 
the solid state electronics (such as increasing of 
performance, reliability and density of elements 
of integrated circuits: diodes, field-effect and 
bipolar transistors) are intensively solving [1- 
6]. To increase the performance of these devices 
it is attracted an interest in the determination of 
materials with higher values of charge carrier’s 
mobility [7-10]. One way to decrease the 
dimensions of the elements of integrated 
circuits is to manufacturing them in thin film 
heterostructures [3-5,11]. In this case, it is 
possible to use inhomogeneity of 
heterostructure and necessary optimization of 
doping of electronic materials [12,13] and 
development in epitaxial technology to improve 
these materials (including analysis of mismatch 
induced stress) [14-16]. An _ alternative 
approaches to increase dimensions of 
integrated circuits are using of laser and 
microwave types of annealing [17-19]. 


Framework, in this article, we introduce an 
approach to manufacture of field-effect 
transistors. The approach gives a possibility to 
decrease their dimensions with increasing their 
density framework, a wideband current 
amplifier with DC-offset cancellation. We also 


consider the possibility to decrease the 
mismatch-induced stress to decrease the 
quantity of defects generated due to the stress. 
We also consider a heterostructure which 
consist of a substrate and an epitaxial layer (see 
Fig. 1). We consider a buffer layer between the 
substrate and the epitaxial layer. The epitaxial 
layer includes into itself several sections, which 
were manufactured by using other materials. 
These sections have been doped by diffusion or 
ion implantation to manufacture the required 
types of conductivity (p or n). These areas 
became sources, drains and gates (see Fig. 1). 
After doping, it is required to annealing of 
dopant and/or radiation defects. The main aim 
of the present paper is to analysis the 
redistribution of dopant and radiation defects to 
determine conditions, which correspond to 
decreasing of elements of the considered filter 
and at the same time to increase their density. 
At the same time we consider a possibility to 
decrease mismatch-induced stress. 


Method of solution 
To solve our aim we determine and analyzed 
the spatio-temporal distribution of 


concentration of dopant in the considered 
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the following form [1,20-24] 
O OC\(x, y, z,t 
poclsy,zt) 


heterostructure. We also determine the 
distribution by solving the second Fick's law in 


OC(x,y,z,t)_ O poClsy.2t) 6) poClay.2t) 
Ot Ox Ox Oy Oy Oz Oz 
Pree sy Lu (x,y, 20) (C(x, y,W,t)dW |+ 
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Oy kT 0 
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Fig. 1a Structure of the considered amplifier [15] 


Substrate 


Fig. 1b Heterostructure with a substrate, epitaxial layers and buffer layer (view from side) 


Here C(xy,z,t) is the spatio-temporal L, 
distribution of concentration of dopant; Q is the surficial gradient; Ic (x, V2.1 )d Z is the 
0 


atomic volume of dopant; Vs; is the symbol of 
surficial concentration of dopant on interface 
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between layers of heterostructure (in this 
situation we assume, that Z-axis is 
perpendicular to interface between layers of 
heterostructure); u4(x,y,z,t) is the chemical 
potential due to the presence of mismatch- 
induced stress; D and Ds are the coefficients of 
volumetric and surficial diffusions. Values of 
dopant diffusions coefficients depends on 


Y 
D, = D,(x, pat) tog SE HE 14 Vemzt) , 


Pp’ ee y, z,T) 


CA 2 
Ds =D, (x, ys aig aes 


Here, D, (x,y,z,T) and Dis (x,,z,T) are the 
spatial (due to accounting all layers of 
heterostruicture) and temperature (due to 
Arrhenius law) dependences of dopant diffusion 
coefficients; T is the temperature of annealing; P 
(xy,z,T) is the limit of solubility of dopant; 
parameter y depends on properties of materials 
and could be integer in the following interval y 
€[1,3] [22]; V (x,y,z,t) is the spatio-temporal 
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properties of materials of heterostructure, 
speed of heating and cooling of materials during 
annealing and spatio-temporal distribution of 
concentration of dopant. Dependences of 
dopant diffusions coefficients on parameters 
could be approximated by the following 
relations [22-24] 


V" (x, y, Zt) 
* 2 (v'y , 
V (x, y, z, V(x, y, Z 
| 1+¢, BBO) aa ate 


distribution of concentration of radiation 
vacancies; V’ is the equilibrium distribution of 
vacancies. Concentrational dependence of 
dopant diffusion coefficient has been described 
in details in [22]. Spatio-temporal distributions 
of concentration of point radiation defects have 
been determined by solving the following 
system of equations [20,23,24] 
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with boundary and initial conditions 
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OT x; 9,0) _¢ OI(x, y,z,t) _0 OI(x, y,Z,t) -0 
Ox a Ox sone Oy ~ 
OI (x, y,z,t) =, 
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OZ A 
I (x,y,Z,0)=fi (x,y,z), V (xy,2,0)=fv (x,),Z). (4) 
Here I! (x,y,z,t) is the spatio-temporal this book); kiv(x,y,z,T), kii(xy,z,T) and 


distribution of concentration of radiation 
interstitials; /* is the equilibrium distribution of 
interstitials; D)(x,y,z,T), Dv(xy,Z,T), Dis(xy, Z,T), 
Dys(x,y,Z,T) are the coefficients of volumetric 
and surficial diffusions of interstitials and 
vacancies, respectively; terms V2(x,,z,t) and 
[2(xy,z,t) correspond to generation of 
divacancies and diinterstitials, respectively (see, 
for example, [24] and appropriate references in 


kyv(xy,z,T) are the parameters of recombination 
of point radiation defects and generation of 
their complexes. 


Spatio-temporal distributions of divacancies 
@y (x,y,z,t) and diinterstitials ®; (x,y,z, t) could 
be determined by solving the following system 
of equations [20,23,24] 
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D : 
+2] dg (yar) E250 E evs Varta fs 
z 0 


Oz ax| kT 
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D (x,y,2,0)=far (x,y,Z), Ov (x,y,Z,0)=fov (x,y,2Z). 


Here Do/(x,y,Z,T), Dov(Xy,Z,T), Dais (X,y,Z,T) and Dovs(x,y,z,T) are the coefficients of volumetric and 
surficial diffusions of complexes of radiation defects; ki(x,y,z,T) and ky(x,y,z,T) are the parameters of 
decay of complexes of radiation defects. 


Chemical potential z in Eq.(1) could be determine by the following relation [20] 
M=E(Z)Qo; [ui(xy,Z,0)+uji(xy,Z,0)]/2, (7) 
1 O Uu,; O Uj 
+ 
2\ Ox, Ox, 
tensor; uj, uj; are the components u;(x,y,z,t), uy(xy,z,t) and u,(x,,z,t) of the displacement vector 
u (x, y, Z, t); X;, Xj are the coordinate x, y, z. The Eq. (3) could be transform to the following form 


where E(z) is the Young modulus, oj is the stress tensor; U = is the deformation 
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where o is Poisson coefficient; a = (as- 


aet)/ae, is the mismatch parameter; as, az, are 
lattice distances of the substrate and the 


36 |-K@)A Ts y.21)-T] 4, SEC) 


expansion; T; is the equilibrium temperature, 
which coincide (for our case) with room 
temperature. Components of displacement 


epitaxial layer; K is the modulus of uniform 
compression; £ is the coefficient of thermal 


vector could be obtained by solution of the 
following equations [25] 
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2 
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heterostructure, 6; Is the Kronecker symbol. With account the relation for oj last system of equation 
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Conditions for the system of Eq. (8) could be written in the form 
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We determine spatio-temporal distributions function corrections [26]. Previously we 
of concentrations of dopant and radiati-on transform the Eqs.(1), (3) and (5) to the 
defects by solving the Eqs.(1), (3) and (5) following form with account initial distributions 
framework standard method of averaging of of the considered concentrations 
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Further, we replace concentrations of dopant values ap. In this situation we obtain equations 
and radiation defects in right sides of Eqs. (1a), for the first-order approximations of the 
(3a) and (5a), on their not yet known average required concentrations in the following form 
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Integration of the left and right sides of the Eqs. (1b), (3b) and (5b) on time gives us possibility to 
obtain relations for above approximation in the final form 
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We determine average values of the first-order approximations of concentrations of dopant and 
radiation defects by the following standard relation [26] 
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We determine approximations of the second 
and higher orders of concentrations of dopant 
and radiation defects framework standard 
iterative procedure of method of averaging of 
function corrections [26]. Framework this 
procedure to determine approximations of the 
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n-th order of concentrations of dopant and 
radiation defects we replace the required 
concentrations in the Eqs. (1c), (3c), (5c) on the 
following sum Qnptp ni (X,y,Z,t). The 
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transformation of the appropriate equations 
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Further we determine solutions of Eqs.(8), i.e. 
components of displacement vector. To 
determine the first-order approximations of the 
considered components framework method of 
averaging of function corrections we replace the 
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Framework this paper we determine the 
concentration of the dopant, concentrations of 
radiation defects and components’ of 
displacement vector by using the second-order 
approximation framework method of averaging 
of function corrections. This approximation is 
usually enough good approximation to make 
qualitative analysis and to obtain some 


quantitative results. All the obtained results 
have been checked by comparison with results 
of numerical simulations. 


DISCUSSION 

In this section, we analyzed dynamics of 
redistributions of dopant and radiation defects 
during annealing and under influence of 


Pankratov., 


International Journal of Advanced Science and Engineering 


www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.9 No.4 3047-3071 (2023) 3066 


mismatch-induced stress. Typical distributions 
of concentrations of dopant in heterostructures 
are presented on Figs. 2 and 3 for diffusion and 
ion types of doping, respectively. These 
distributions have been calculated for the case, 
when value of dopant diffusion coefficient in the 
epitaxial layer is larger, than in the substrate. 
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The figures show, that inhomogeneity of 
heterostructure gives us possibility to increase 
compactness of transistors. At the same time 
one can find increasing homogeneity of dopant 
distribution in doped part of epitaxial layer. In- 
creasing of compactness of transistors gives us 
possibility to increase their density. 
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Fig.2: Distributions of concentration of infused dopant in heterostructure from Fig. 1 in direction, 
which is perpendicular to interface between epitaxial layer substrate. Increasing of number of curve 
corresponds to increasing of difference between values of dopant diffusion coefficient in layers of 
heterostructure under condition, when value of dopant diffusion coefficient in epitaxial layer is larger, 
than value of dopant diffusion coefficient in substrate 
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Fig.3: Distributions of concentration of implanted dopant in heterostructure from Fig. 1 in direction, 
which is perpendicular to interface between epitaxial layer substrate. Curves 1 and 3 corresponds to 
annealing time © = 0.0048(L,2+L,2+L,2)/Do. Curves 2 and 4 corresponds to annealing time © = 
0.0057(L,2+Ly2+L72)/Do. Curves 1 and 2 correspond to homogenous sample. Curves 3 and 4 
corresponds to heterostructure under condition, when value of dopant diffusion coefficient in 
epitaxial layer is larger, than value of dopant diffusion coefficient in substrate 
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Fig. 4: Spatial distributions of dopant in heterostructure after dopant infusion. Curve 1 is idealized 
distribution of dopant. Curves 2-4 are real distributions of dopant for different values of annealing 
time. Increasing of number of curve corresponds to increasing of annealing time 
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Fig. 5: Spatial distributions of dopant in heterostructure after ion implantation. Curve 1 is idealized 
distribution of dopant. Curves 2-4 are real distributions of dopant for different values of annealing 
time. Increasing of number of curve corresponds to increasing of annealing time 


The second effect leads to decreasing local 
heating of materials during functioning of 
transistors or decreasing of their dimensions for 
fixed maximal value of local overheat. However 
framework this approach of manufacturing of 
bipolar transistor it is necessary to optimize 
annealing of dopant and/or radiation defects. 
Reason of this optimization is following. If 
annealing time is small, the dopant did not 
achieve any interfaces between materials of 
heterostructure. In this situation one cannot 


L,LyL, 


find any modifications of distribution of 
concentration of dopant. If annealing time is 
large, distribution of concentration of dopant is 
too homogenous. We optimize annealing time 
framework recently introduces approach 
[15,25-32]. Framework this criterion we 
approximate real distribution of concentration 
of dopant by step-wise function (see Figs. 4 and 
5). Farther we determine optimal values of 
annealing time by minimization of the following 
mean-squared error 


—— fjffe (x, y,z,0)-y (x, y,z)|dzd ydx, (15) 
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where y (x,y,z) is the approximation function. 
Dependences of optimal values of annealing 
time on parameters are presented on Figs. 6 and 
7 for diffusion and ion types of doping, 
respectively. It should be noted, that it is 


necessary to anneal radiation defects after ion 
implantation. One could find spreading of 
concentration of distribution of dopant during 
this annealing. In the ideal case distribution of 
dopant achieves appropriate interfaces between 
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materials of heterostructure during annealing of dopant. In this situation optimal value of 
radiation defects. If dopant did not achieve any additional annealing time of implanted dopant 
interfaces during annealing of radiation defects, is smaller, than annealing time of infused 
it is practicably to additionally anneal the dopant. 
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Fig.6: Dependences of dimensionless optimal annealing time for doping by diffusion, which have been 
obtained by minimization of mean-squared error, on several parameters. Curve 1 is the dependence of 
dimensionless optimal annealing time on the relation a/L and €= y= 0 for equal to each other values 
of dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of 
dimensionless optimal annealing time on value of parameter ¢for a/L=1/2 and €= y= 0. Curve 3 is the 
dependence of dimensionless optimal annealing time on value of parameter € for a/L=1/2 and e= y= 
0. Curve 4 is the dependence of dimensionless optimal annealing time on value of parameter y for 
a/L=1/2 and e= €=0 
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Fig.7: Dependences of dimensionless optimal annealing time for doping by ion implantation, which 
have been obtained by minimization of mean-squared error, on several parameters. Curve 1 is the 
dependence of dimensionless optimal annealing time on the relation a/L and €= y= 0 for equal to each 
other values of dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence 
of dimensionless optimal annealing time on value of parameter ¢ for a/L=1/2 and €= y= 0. Curve 3 is 
the dependence of dimensionless optimal annealing time on value of parameter & for a/L=1/2 and ¢= 
y= 0. Curve 4 is the dependence of dimensionless optimal annealing time on value of parameter y for 
a/L=1/2 and ¢= €=0 
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Fig. 9: Normalized distributions of charge carrier mobility in the considered heterostructure. Curve 1 
corresponds to the heterostructure, which has been considered in Fig. 1. Curve 2 correspond to a 
homogenous material with averaged parameters of heterostructure from Fig. 1 


Further, we have analyzed the influence of 
relaxation of mechanical stress on distribution 
of dopant in doped areas of heterostructure. 
Under the following condition #< 0 one can find 
compression of distribution of concentration of 
dopant near interface between materials of 
heterostructure. Contrary (at 4>0) one can find 
spreading of distribution of concentration of 
dopant in this area. This changing of 
distribution of concentration of dopant could be 
at least partially compensated by using laser 
annealing [29]. This type of annealing gives us 
the possibility to accelerate diffusion of dopant 
and another process in the annealed area due to 
inhomogenous distribution of temperature and 
Arrhenius law. Accounting relaxation of 


mismatch-induced stress in heterostructure 
could leads to changing of optimal values of 
annealing time. Mismatch-induced stress could 
be used to increase density of elements of 
integrated circuits. On the other hand could 
leads to generation dislocations of the 
discrepancy. Fig. 8 shows distributions of 
component of displacement vector, which is 
perpendicular to interface between layers of 
heterostructure. 


CONCLUSION 

In this paper, we modeled the redistribution of 
infused and implanted dopants with account 
relaxation mismatch-induced stress during 
manufacturing field-effect transistors 
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framework a wideband current amplifier with 
DC-offset cancellation. We formulate 
recommendations for optimization of annealing 
to decrease dimensions of transistors and to 
increase their density. We also formulate the 
recommendations to decrease mismatch- 
induced stress. Analytical approach to model 
diffusion and ion types of doping with account 
concurrent changing of parameters in space and 
time has been introduced. At the same time the 
approach gives us possibility to take into 
account nonlinearity of considered processes. 
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